Translational invariance of Coulomb series and symmetric potentials in crystals 
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It is shown that Coulomb series are to be considered within a special mode of summation so as to 
describe bulk properties of crystals. The translational invariance is then an explicit integral property 
of Coulomb series that is tantamount to the effect of invariant periodic boundary conditions discussed 
earlier. Absolute bulk potentials with zero mean value are then substantiated as a unique solution 
in the general case of triclinic lattices. An invariant treatment of the bulk Coulomb energy follows 
therefrom. The potential symmetry is verified for simple point-charge lattices and is connected with 
the centre of gravity of the potential field that is relevant to non-local charges as well. 
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I. INTRODUCTION 



Sums of Coulomb potentials over lattices (lj are basic 
for describing the cohesive energy of solid state g, [| ||, 
H ^, 0, §]■ Such sums also determine local potentials 
responsible for the electronegativity and charge transfer 
in crystals §, |, [l0[ || HI HE HI and so specify 




36[ |7|, H H, li- 



the electronic band structure 
[H],f2| H H |§ [§,[27j. 

The problem which still exists is to determine bulk 
Coulomb series so that the result of summation in all the 
above events be unique in every particular case. The po- 
tential symmetry is then evident at least in the simplest 
point-charge lattices of AB type [g || g {$G|, [H], 
Therefore, the bulk Coulomb energy can be specified by 
particular values of the absolute local potentials therein 

n m in" 



How- 
ever, different modes of summation with various types 
of parametrization either in the direct space p8| , |3l| , 34], 
H H, H |J, HI Pl ^ M |50| o r with making use 
of the reciprocal one |2j, |33|, |38|, g |5l|, ||, || |§ || Q 
lead to controversial results. This fact is due to the condi- 
tional character of convergence of infinite Coulomb series 

To be unique, the bulk solution of interest must, cer- 
tainly, be independent of the optional parameters of the 
mode of summation ||| [|| Q ||, As a hint 

about it, surface effects are to be irrelevant to bulk prop- 
erties §, |J ||, || H g ||, £§ || g ^, ||. 
All the statistical treatment of solids is based on this 
statement ||| U, Hi- Zero mean bulk potential is 
declared therefrom regardless of the summation scheme 
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^|,£0^, |3J, |38|, ^3|, ^6, 

the case of parent infinite series 
periodic boundary conditions 



But it is neither 
381 |69[] nor the case of 
|7Q ] while they may be 
modified by a constant potential as an available periodic 
solution of Laplace's equation jl8|, |2t| [43|, ^] . 

As found jn] f72|, the solution of interest arises under 
periodic boundary conditions invariant to the definition 
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of the unit cell relative to the crystal structure, which 
generalize the translational invariance based on periodic 
replications of a certain basic unit cell j|, ||(], |||, 72, 
[74], f75|. This result substantiates Ewald's proposal 1 32 ] 
to exclude the uniform potential component [^9[ 74] 
as well as the energy effect of surface polarization 72, 
[74] . [tg| , [77t from the results of summation that was not 
convincing for a long time § [l§ |fj, ^ [7§ |Tj, 



Moreover, though the energy value for a neutral unit 
cell in the bulk is formally independent of any conceiv- 
able change in the potential origin || [H], p2l [fj], |4j , 
the bulk Coulomb energy and absolute bulk potentials 
appear to be functionally connected. The fact that elec- 
trostatic quantities in the bulk are independent of real 
surfaces implies that the configuration of any surface is 
to be consistent with the bulk state || ^, ^2|, |6l], [7(], 
pE |B6| |87| , |8S| , p9[ . In particular it means that the surface 
conditions have to reproduce the bulk state in question 



It is important that the same problem is relevant to 
the definition of direct bulk Coulomb series as such. In- 
deed, pristine Coulomb series, albeit expressed in terms 
of periodic charge distributions, do not contain the peri- 
odicity as a property of summation in an explicit fashion 
p3| . In the present paper we deduce that invariant peri- 
odic boundary conditions actually imply a special mode 
of summation as a continuous procedure with piecewise 
smooth intermediate boundary surface and with the fur- 
ther averaging of the oscillating result over the period 
specified by lattice translations. As shown, apart from 
the convergence conditions, an additional one should be 
imposed on the charge distribution per unit cell so as to 
represent Coulomb series in a conventional form. As an 
example, the symmetric potentials in cubic crystals are 
readily calculated in an effective manner. The topologi- 
cal nature of the potential symmetry is discussed. 



2 



II. AMBIGUITY ABOUT ABSOLUTELY 
CONVERGENT COULOMB SERIES 

If p(r) is a certain charge distribution attributed to 
a unit cell, then the direct Coulomb sum associated 
with p(r) and describing the electrostatic potential at 
the point r can be written as 



U cd (r) = 



p(r') dr' 
-r' R, 



(1) 



where i runs over the Bravais lattice composed of the 
unit-cell origins R;, the integral is over the volume V 
occupied by p(r'), which can spread beyond the unit-cell 

volume §, [21], |2§ §§, ^ §1 13 ' the P rime on the sum " 
mation sign implies missing the contribution of any point 
charge if it happens at r. For the absolute convergence 
of ( |l|), p(r) is supposed to obey the following constraints 
|Jf3 [4§ g g ||, g ^, 0, 0: 



/ p(r) dr = 0, 

M„ ee / r M p(r) dr = 0, 
Jv 

G/j,v = / r^r v p{r) dr = at ^ ^ ^, 



G xx — Gyy — G zz — H, 



(2) 

(3) 

(4) 
(5) 



where are Cartesian co-ordinates of r, H is a certain 
constant. Condition (|[) ensuring the local neutrality is 
conventional. Relations (HWpf ) can be maintained by 
symmetry |, H pi, |a |T |1. If they do not hold 
for any initial p 1IU (r) [§ §§f g El' 
p lm (r) can always be modified properly ||7L f72||~^ The 
implication is that the electrical neutrality is the only 
fundamental constraint in crystals || |2(], |33[ [35l |5f], 
p^ . [7G|, |80| , |Sl[ |93[ ] , but a rather arbitrary charge distribu- 
tion p m ^r) transformed into p(r) driven by (||)-(||) with 
an optional H may be adopted for the purpose of the 
lattice summation specified as a definite procedure in (|l|) 

The residual effect of H is evident from the mean po- 
tential value Ucd associated with (Q) and it is indepen- 
dent of the order of summation due to (Q)-(|5|). Following 
Bethe |94|, we have 



U C d = 



2tt / 2 , 2irH 

6V I V 



(6) 



where v is the volume of the unit cell, the left equality of 
(§) is typical §, |(| |§ frl || |J, |9| and is further 
converted with the help of (jjf). We see that the value of 
Ucd turns out to be optional through H. This inference 
can be naturally extended to the potential C/cd(r) as such 
po| , |69j . However, it is also clear that the effect of H is 
artificial and should be removed. 

Although the case of H = is of great importance, as 
will be clear later on, it does not exhaust the situation 



at hand. To gain insight into this problem, we turn to 
the limiting behaviour of (|l|) upon the direct summation 
over closed shells of one unit cell thickness, enveloping 
a reference unit cell in a consecutive fashion. According 
to |§)-(§), such an order of summation is not more de- 
manded by the arguments of convergence [pp|, pq , p8[ |35| , 



51J, |52|, £5[|8|, |60|, [73|, |7§ |9(f, but al- 
lows for the equal uniformity along each crystallographic 
axis 

@ H H H !§• The absolute convergence of 
(0) implies that the contribution of remote shells with 



numbers m > m c 3> 1 is negligible: 



lim 



p{r) dr 
IRi_ +r 



0. 



(7) 



where i m runs over the unit cells belonging to the mth 
shell. Relation (Q) persists as long as remote unit 
cells arc added as permanent bricks described by p(r) 

|, u n 10 m M> D> H H H H H 



Since 

individual features of such bricks have to be lost in the 
crystal interior [^0[ |44|, |59|, the fortuitous contribu- 
tion of H to the lattice sum is patently generated by 



outer shells of summation 
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and so must be 



eliminated as irrelevant to bulk properties now. 



III. CONTINUITY OF SUMMATION AND 
INVARIANT PERIODICITY 

Inasmuch as (^) holds due to the discrete character of 
summation over shells there, it is reasonable to modify 
this type of summation so as to remove that discrete- 
ness, but without destruction of the charges taken into 
account ||(], |it], |(3(| |93j. To this end, we conceive that 
our crystal is enclosed in an imaginary large bounding 
box homothetic to the unit cell, with a reference unit 
cell at its centre, and the summation region is restricted 
by the box volume. The summation over shells is then 
treated as a continuous procedure of enlarging this box. 
In a triclinic lattice with the lattice parameters a, b and 
c along unit vectors e a , e;, and e c , the evolution of the 
bounding box can be described by the parameters A, B 
and C corresponding to e a , and e c and increasing lin- 
early in an evolution parameter / at 



ABC 
a b c 



(8) 



as shown schematically in Fig. |l|(a). 

The event of p(r) arranged within the unit-cell par- 
allelepiped is first considered. The limiting procedure in 
(0) can then be readily associated with the bounding-box 
evolution if we imagine that the m c th shell is truncated 
by the bounding box at a given moment. It implies that 
the inner part of this shell gives an additive contribution 
to the maternal sum over m < m c and this shell is re- 
moved from (Q) and replaced by the next one when its 
outer boundary is achieved by the bounding box. The 
periodicity of the latter events upon the evolution in 
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*c/)-*] E 

;e{A} 



_LA 



where t(/) = af — a/2, r = te a +£>e& + ue c , 

i? = [l — cos 2 a — cos 2 (3 — cos 2 7 

1/2 

+ 2 cos a cos /? cos 7] 



(9) 



(10) 



at (e a et) = cos a, (e;,e c ) = cos/3 and (e c e a ) = cos 7. 

The averaging of $a(/) over the period of oscillations 
is readily defined in the form 



FIG. 1: Summation over the shells as a uniform continuous 
procedure: (a) Remote unit cells, small regular parallelo- 
grams, are united into the shells bordered by heavier lines, 
with the hatched area of an incomplete outermost shell trun- 
cated by the bounding box depicted by the heaviest lines. The 
bounding-box interior contributing to the lattice sum enlarges 
uniformly upon increasing the evolution parameter /. (b) Os- 
cillations of the limiting value of a certain inner potential U 
as a function of / are expected relative to its magnitude spec- 
ified by complete outermost shells. This initial magnitude 
marked by the filled circles is treated as an origin, for conve- 
nience. Averaging over oscillation periods shown by vertical 
dotted lines converts the oscillating U into an invariant value 
depicted by the heavy horizontal line, of which position de- 
scribes the potential shift of interest. 



question reflects the translational invariance of (0) ac- 
companied by some periodic change in the value of the 
resulting maternal sum, as sketched in Fig. 0(b). This 
limiting effect can naturally be analysed within a single 
period of oscillations at < / < 1, so that / becomes a 
measure of the incorporated share of the m c th shell. 

Note that now the expected result depends on our sum- 
mation mode exhibiting the translational invariance as 
an integral property of Coulomb series |33[ ]. To make 
the above oscillations surely finite and minimal, all the 
faces and edges of the bounding box are assumed to be so 
charged as to retrieve local constraints (0) and (0) for the 
truncated unit cells at every moment of the bounding- 
box evolution. Then the potential effect of the remote 
charged edges of the truncated m c th shell modified by 
charges on the corresponding edges of the bounding box 
vanishes because it is specified by the mean charge den- 
|, which are zero by definition, 
ration of the modified m c th shell 
faces to (0), we focus on the Ath face disposed in the e a 
direction, with R^ describing the unit cells on this face. 
Here we are interested in the potential contribution at 
the centre of the bounding box. In the expansion of the 
integrand in (0) in powers of r/Ri, the ratio of moduli, 
only terms linear in the components R^~ A of Ri along an 
outward normal to the Ath face survive |(il], [si, . 
Their contribution to (m) can be written as 



sities [|28|, |59j |60] |6J |8£ 

To discuss the contri 



Q2 rt(f) rb/2 rc/2 

$A(f) = ^—3 I dt dp dup(t,p,u) 

Smfl J_ a / 2 J-b/2 J-c/2 



*a(/)# 



(11) 



and is the final logical step of elimination of the summa- 
tion boundary, because the boundary which has neither 
definite structure nor definite position does not more ex- 
ist. This case is naturally associated with the topologi- 
cally endless lattice. Indeed, the same result is valid for 
the contribution <?_a of the —Ath face of the m c th shell 
by symmetry, but their combination may be represented 
as 



1 . 



*a(/) + *-a(1-/) 



df = 2$ A , 



(12) 



where the integrand corresponds to the closure of the 
lattice in the e a direction along the ±Ath bounding-box 
faces, in which case these auxiliary bounding-box faces 
merge and so neutralize each other by definition, but pe- 
riodic boundary conditions arise instead. The averaging 
over / in (|l^) then implies the invariance of the periodic 
boundary conditions to the optional position of the plane 
of joining inside particular unit cells |0, |7/| . 

On determining <^b and <Pc in the same manner and 
taking our previous result j7l], [72| into account, for the 
total potential eliminating the boundary effect we obtain 



ftop 



2($, 



2irH 



(13) 



Note that issue (|T^) is relevant to the foregoing periodic 
boundary conditions, which hold for each basic crystal- 
lographic direction and so are three-dimensional. The 
duality of the treatment of the procedure of averaging 
makes this inference geometrically noncontradictory. 

It is significant that like (0) , result ( |l3"| ) is independent 
of the particular features of the crystal symmetry. Ow- 
ing to the obvious additivity of the potential effect, result 
(O) can be extended to overlapping p(r), as well as to 
more complicated Bravais lattices of higher symmetry. 
By induction, the validity of (|l^) can also be extended 
to events of p(r) spreading to infinity as long as the ba- 
sic statements that H is finite and that the stationary 
limiting regime is provided by (0) are maintained. As a 
result, the general concept that a rather arbitrary ^'"'(r) 
can be involved for describing a crystal ^0), ^8) is 

substantiated. 
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TABLE I: Tentative unit-cell point charges Sj in units of k for each of til primitive cubic lattices, with the charge positions 
generated from those indicated, in units of a, in the parentheses by virtue of cubic transformations. 



Crystal 


n L 


si (0,0,0) 


s 2 ( 4 j 4 j 4 ) 


s 3 (i,0,0) 


S4 (i,i,0) 


s 5( §1414) 


s 6( 2 , 2 , 2 ) 


S7(l,0,0) 


„ (3 3 3\ 


1 3 3 3 \ 


CsCl 


1 


337 
256 








— 





513 
4096 


27 
512 





1 

4096 


NaCl 


4 


32 
23 




7 

4C 


3 

92 




1 

92 








ZnS 


4 


415 
384 


3 a 
IB 




1 
144 


1 a 
48 




1 

2304 






Cu 2 


4 6 


17 

8 


3 c 
8 




3 

32 


3 c 
64 






1 c 
64 




CaF 2 


4 


35 
16 


49 
256 




1 

64 


5 
256 










BiF 3 


4 


591 
184 


1 
4 


61 
368 


13 
736 




1 

1472 









"These charges are at the proper tetrahedral positions only. 

'The final structure is made up of primitive lattices rotated about a fourfold symmetry axis. 
c These charges are on (1, 1, 1) diagonals only. 



TABLE II: As a result of the direct summation over point-charge cubic lattices, the absolute bulk potentials Ubj, in units of 
ft/a,at the unit-cell symmetric points specified in the parentheses and the Madelung constant q:m (the bulk energy £ b described 
by (|16[) in units of —k 2 /a). The positive ion is at the (0, 0, 0) point. The repeating potential values are indicated in symbols. 



Crystal 


U bl (0,0,0) 


Ub2 (j,J)f) 


Uh3 (—4,474) 


(7b4(i,0,0) 


u b5 (5,1,0) 


Ub6 ( 2 ' 2 ' 2 ) 




CsCl 


-2.03536151 


0" 


u b2 a 


0.48658923 a 


-0.48658923 a 


2.03536151 


2.03536151 


NaCl 


-3.49512919 


0" 


u b2 a 


3.49512919 


t/bi 


u b4 


3.49512919 


ZnS 


-3.78292610 


3.78292610 


0.28779691 " 


-0.28779691 a 


t/bi 


u b4 a 


3.78292610 


Cu 2 


-3.78292610 


6.47653093 


0.55497170 a 


-0.28779691 a 


u bl 


u bA a 


10.25945703 


CaF 2 


-7.56585221 


4.07072302 


u b2 


-0.57559383 a 


u bl 


u b4 a 


11.63657523 


BiF 3 


-11.06098140 


4.07072302 


u b2 


2.91953536 


u bl 


u b4 


22.12196279 



"The corresponding lattice positions are interstitial. 



Note that our procedure of averaging the potential ef- 
fect of the bounding-box evolution resembles antecedent 
approaches @ |]J || |f| || || |§, which might be 
associated with the averaging over a symmetric set of 
discrete points on the oscillatory potential curve. 

According to ([!]) and ([ll]), the bulk potential field in a 
crystal is as follows: 



U h (r) = U cd (r) + 



top ■ 



(14) 



Although relationship ( Jl4|) is anticipated J3^, pSft , the 
topological treatment of $top 0, Q differs from earlier 
interpretations (l9[ 0. On combining ( [Ti]) and (^|), for 
the mean bulk potential we obtain the expected result 
Ub = 0. Hence, Ub(r) defined by (jTJ) is independent of 
H and does not contain any uniform component. The 
potential symmetry follows therefrom. 



IV. DISCUSSION 



We may apply ( |l3| ) and (14) to the rapidly conver- 
gent direct Coulomb summation for several cubic 
point-charge structures with the lattice spacing a and 
with charges measured in units of k, the modulus of a 
minimal point charge. These structures are treated as 



decomposed into n L primitive lattices (|, Q with the 
unit-cell charges compiled in Table Q. According to (|lJ 
for each structure in question we obtain 



TTK 

^top = X { 

a 



'4(s 6 


f s 7 + 


9s 9 ) 




CsCl , 


4(s 3 


f 4s 4 


+ 456) 




NaCl , 


2(s 2 


4- 8s 4 


+- H«6 


+ 8s 7 ) 


ZnS , 


(S2 + 


I6S4 


4-llSS 


+ 9s 8 ) 


Cu 2 0, 


4(s 2 


f 4s 4 


f lls 5 


) 


CaF 2 , 


A(s 2 


f S3 + 


4s 4 + 


4s 6 ) 


BiF 3 



(15) 



as the summary topological contribution of all til sublat- 
tices. The resulting numerical values calculated through 
(0) for each local potential specified by ( 3) at some sym- 
metric points are arranged in Table p , wh ere ( |l5| ) re- 
trieves the potential symmetry |l|, ||, |3C , |lOCf| . 

Of course, dealing with an arbitrary p lnl (r), it is 
reasonable to go over to p(r) specified by H = 0. 
Then the boundary effects are formally excluded and 
the analysis is available in terms of direct sums only. 
In this case Uh(r) — J7cd(r) that explains Evjen's ap- 
proach j34|] to NaCl, for example. Moreover, the bulk 
Coulomb energy Eh per unit cell can then be confi- 
dently written in terms of ([!]) in a conventional manner 
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1 f 7 , ,^'f p(v')dr' 



r - r' - R, 



= ~Jdrp(r)U h (r) = ±J^dvp ini (r)U b (r), (16) 

in agreement with ([li]) . The last relation in ([l6|) is due to 
the absolute character of the potential field J7b( r ) and, 
in general, implies that every charge distribution con- 
sistent with the crystal structure may be incorporated 
there, if we are interested in the value of £h at given 
?7b(r)- In particular, the Madelung constants in Table |H| 
are obtained with making use of the potentials U^j (r) on 
integer point charges, in agreement with t he c o rrcspond - 
ing precise data M ||, || ||, [lM] |0§ [l0§ pj . 
Bearing in mi nd J^ jT p2[ that the customary functional 
relation §, [l|, 



<5p ini (r) 



f/b(r) 



(17) 



as a base for the self-consistent determi- 
(r)J| Jj, ||, §,§§, any difficulties 
3, |80|, p3[ in defining £^ are obviated via 



holds for (jO 
nation of p" 
g^, g |_ 

The topological effect of <P tov is evident in the hy- 
pothetical case of lattices composed of non-interacting 
neutral atoms, which are to b e spherical and non- 
overlapping for this purpose [105 . In the simplest event 

~ Its' 



0, 52 



106, 107], we consider such atoms as solid 



spheres of radii tq, charged uniformly, with neutralizing 
positive nuclear charges Z at their centres. One atom 
per unit cell is assumed. Then relations (||)-(^) hold by 
symmetry. According to @ and (13), we derive 



<?top = 



27rZr, 



o 



5 1' 



(18) 



On the other hand, the sum in (Fy) is reduced to a sol e 



term if r is inside an atom and is zero otherwise [ 105 
The calculation of the mean value Ucd is then straight- 
forward and we get 



3Z 



(19) 



where U is the mean intratomic potential, the factor 
in the parentheses renormalizes U per v. According to 
@, @ and @, we reach U h = again. However, 



as a centre of gravity [M 13, [L9L K37], <P t op now con 



nects intratomic potentials which dwell as formally in- 
dependent in the case at hand. It is worth noting that 



this effect of ( [181) is absent in low-dimensional lattices 
po[ H?], [|(], [8^, |91] where $ t o P = 0. But such a negative 
shift due to derealization of electrons in atoms and ions 
is a general effect which, apart from the effects of ionic 
overlap [^0[ ^l], |2{| , modifies the three-dimensional point- 
charge lattice potentials. According to estimate (|l8|), this 
shift is compared to potentials in Table || and may ac- 
count for trapping small atoms and ions at interstices 
109, 110, 111] as well as fo r the asymmet ry of 

ToT 



112, 113, 114, 115, 116| 



vacancy concentrations 

Note that every artificial derealization of charges 




f§ p, ||, |?J always pro- 

duces a certain potential shift of the discussed nature, 
with including uniform charge distributions J^, |33|, |6T 
For instance, according to 



33|, |91[. For instance, according to (^J and 
( |13| ) , in the case of simple cubic lattices a uniform charge 
backgroung with the total charge Q per unit cell gives rise 
to <?top = 7rQ/6a. This term compensates the ordinary 
background effect |82[] . However, such contributions can- 
cel out in the superposition if, for example, two identical 
sublattices composed of unlike charges, but with common 
parameters of charge derealization are studied Q. 



V. CONCLUSION 

In summary, a special mode of summation associated 
with invariant periodic boundary conditions must be in- 
cluded in the definition of Coulomb lattice series so as to 
describe bulk properties of crystals. As a result, abso- 
lute electrostatic potentials arise as a unique solution of 
interest. Zero mean bulk potential ensuring this unique- 
ness within every summation scheme devised to date is 
substantiated rigorously. Based on absolute bulk poten- 
tials, the invariant treatment of the bulk Coulomb energy 
is proposed. The limiting cases estimated highlight the 
potential symmetry as a peculiar feature of the three- 
dimensional potential map. 
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